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Abstract. The spatial self-organization of a Bose-Einstein condensate (BEG) in a high-finesse linear optical 
cavity is discussed. The condensate atoms are laser-driven from the side and scatter photons into the cavity. 
Above a critical pump intensity the homogeneous condensate evolves into a stable pattern bound by the 
cavity field. The transition point is determined analytically from a mean-field theory. We calculate the 
lowest lying Bogoliubov excitations of the coupled BEC-cavity system and the quantum depletion due to 
the atom-field coupling. 

PACS. 03.75.Kk Dynamic properties of condensates; collective and hydrodynamic excitations, superfiuid 
flow - 37.10.Vz Mechanical efl^ects of light on atoms, molecules, and ions 



1 Introduction 

Atoms in a high-Q cavity, even in the dilute gas limit, 
form a unique many-body system due to the long-range 
interaction mediated by the commonly coupled cavity field 
mode. We have recently discovered that a gas of laser- 
driven cold atoms scattering light into a cavity produces 
a phase transition between the homogeneous spatial dis- 
tribution and a regular periodic pattern [T]. The ordered 
phase is, in fact, a 'giant' cubic crystal with a lattice con- 
stant of the optical wavelength, which is bound by the ra- 
diation field scattered by the atoms into the cavity. This 
self-organization has been observed in experiment [2^ , and 
has become a workhorse to study many-body phenom- 
ena in the strong coupling regime [3||4][5|. The realization 
of collective effects with material exhibiting wavelike be- 
haviour is within reach since the first strong coupling ex- 
periments in the ultracold atom temperature regime have 
already been reported [Bjt7,8,,9^. A closely related, many- 
atom system in cavity is the correlated atomic recoil laser 
(CARL), demonstrated recently [10], which has also been 
studied in the ultracold gas limit [TT]. CARL gain with 
BEC in free space has been observed |12j . 

In this paper we study the self-organization phenom- 
enon in the case of coupled matter and radiation fields. 
As a first step in understanding the complex dynamics 
of this system, we neglect quantum statistical effects |JJ5i 
m] and resort to a mean field approximation. We adopt 
the model developed by P. Horak et al. [15] for a dif- 
ferent geometry. In their scheme, the coherent laser field 
was injected directly into the cavity. Then the transition 
between homogeneous and periodically modulated densi- 
ties is smooth as a function of the external pump power. 
By contrast, we consider the case of illuminating directly 



the atoms, which yields a phase transition, i.e., an abrupt 
change of the stationary state of the system at a well de- 
fined threshold pump power. Note the delicate difference 
between stationary state and ground state of the system. 
The condensate does not follow a coherent Hamiltonian 
dynamics because of the coupling to the decaying cavity 
field. The total system is dissipative owing to the irre- 
versible photon loss from the cavity. Thus we have to con- 
sider the steady-state of a driven open system far from 
equilibrium and perturbations around this steady-state. 

The paper is organized as follows. In Sec.|2]we present 
the system and the mean field model based on coupled 
Gross-Pitaevskii and cavity field equations. Next, in Sec. [3] 
we calculate the steady-state and show that it undergoes 
a symmetry-breaking transition at a critical pump power. 
The pumping threshold for self-organization will also be 
discussed in terms of collective excitations in Sec. |4j The 
spectrum of low lying excitations is presented in the full 
range of the pumping strength in this section. The strong 
collective coupling regime of cavity QED is marked by a 
different type of transition where the ordered phase per- 
mits defect atoms in the self-organized lattice. This effect 
is presented in Sec. [5] Section [6] is devoted to studying 
quantum depletion which is a crucial issue concerning the 
stability of the condensate in a cavity and also with re- 
spect to the validity of the model. Finally, we conclude in 
Sec.0 



2 Model 

We consider a pure Bose-Einstein condensate (BEC) in- 
teracting with a single-mode of a high-Q optical cavity. 
The condensate atoms are coherently driven from the side 
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by a laser field with frequency oj, directed perpendicu- 
larly to the cavity axis. The laser is detuned far from 
the atomic transition uja, that is, \Aa\ ^ 7, where 27 is 
the full atomic linewidth at half maximum and the atom- 
pump detuning is A a = uj — u>a- This condition ensures 
that the electronic excitation is extremely low in the con- 
densate atoms, hence the spontaneous photon emission is 
suppressed. At the same time, the laser field is nearly res- 
onant with the cavity mode frequency loc, i.e. \Ac\ ~ k, 
where k is the cavity mode linewidth and the cavity-pump 
detuning is Ac = to — uc- The scattering of laser pho- 
tons into the cavity is thus a qusi-resonant process and is 
significantly enhanced by the strong dipole coupling be- 
tween the atoms and the mode due to the small volume of 
the cavity. This coupling strength is characterized by the 
single-photon Rabi frequency g, which is in the range of 
K. Therefore, although the condensate is hardly excited, 
it can efficiently scatter photons into the cavity. 

For the sake of simplicity, we describe the dynamics in 
one dimension x along the cavity axis. The cavity mode 
function is then simply cos kx. This model can apply e. g. 
to a cigar shaped BEC tightly confined in the transverse 
directions by a strong dipole or magnetic trap, so that the 
transverse size of the condensate w is smaller than the 
waist of the cavity field. The pump laser is assumed to be 
homogeneous along the cavity axis therefore it is described 
by a constant Rabi frequency J7. 

The condensate contains a number of N atoms as- 
sumed to have the same wave function |V'(^))- The cavity 
field is assumed to be in a coherent state described by the 
complex amplitude a. These approximations imply that 
the quantum state of the system is factorized: entangle- 
ment between the condensate and the cavity field [511TB] 
is neglected, which can be done for large enough cavity 
photon number 

The cavity field is subject to the strong refractive in- 
dex effect of the optically dense condensate. At the same 
time, the evolution of the condensate wave function is de- 
scribed by a Gross-Pit aevskii- type equation, including the 
mechanical effect of the radiation field in the cavity. The 
system of coupled mean-field equations is 

i^a - [-Ac + N{U{x)) -iK]a + N^x)), (la) 
ot 



z|^(a:,t) = {2£^ + |a(i)|2[/(x) + 2Re{a(t)K(x) 

+Ng,mx,t)\^}ij{x,t). (lb) 

These equations are closely related to the semiclassical 
model for pointlike atoms [T| fT7l[T8] . and rely on the same 
approximations as the ones used in the cavity-pumping 
geometry [15]. Alternatively, this is the single-mode field 
case of the model of Ref. [19] which relies on a canoni- 
cal Lagrangian approach to the action-backaction in the 
coupled atom- light system in the large detuning limit. 

Let us discuss the physical meaning of the terms cou- 
pling the photon and the matter fields. Each atom shifts 
the cavity resonance frequency in a spatially dependent 



manner by U (x) = Uq cos^ kx. The maximum shift is Uq = 
g'^/AA, obtained at the antinodes of the mode function. 
In the mean-field approximation, the shift has to be spa- 
tially averaged over the single-atom wave function, giving 
the frequency shift per atom {U{x)) = Uq {ipl cos^ kx 
It is worth noting that, even if the one-atom light shift Uq 
is small compared to the cavity decay k, one can achieve 
the interesting strong collective coupling regime of cav- 
ity QED, N\Uo\ > K, with a trapped BEC. Feeding the 
cavity by laser scattering on the atoms appears as an ef- 
fective pump with strength 77 = Qg/AA- This process has 
a spatial dependence inherited from the mode function, 
rit{x) — rj cos kx, and this term also has to be averaged 
over the condensate wave function. 

The well known Gross-Pitaevskii equation (GPE) in 
one dimension in Eq. ( lb) describes the motion of a scalar 



condensate induced by the optical forces in the cavity field. 
The back action of the light shift is the term proportional 
to U{x) = UqCos^ kx and the cavity photon number \a\^. 
This term is periodic with half of the wavelength A and, 
in free space, is referred to as an "optical lattice". The 
back action of the coherent scattering of photons between 
the transverse pump and the cavity mode is the term pro- 
portional to rit{x) = Tycos kx and has a periodicity of A. 
The last term of the GPE accounts for the s-wave col- 
lision of the atoms, its strength is related to the s-wave 
scattering length a by = 47r/ia/(mu'^), and depends on 
the transverse size of the condensate w. 

In this paper we will consider the case Uq < 0, i. e. 
large red detuning, where the atoms behave as high field 
seekers. Consequently, for nonzero cavity field the con- 
densate atoms tend to localize around the field antinodes, 
thereby maximizing their coupling to the light field. 



2.1 Scaling with the atom number 

It can be checked that Eqs. ([l]) with the field amplitude 
replaced by a/^/N are invariant under the scaling of the 
parameters such that NUq, Ng^, and ^fNrj is kept con- 
stant. That is, in the mean field model, the atom number 
can be incorporated in the system parameters and the 
field amplitude variable (with the proposed rescaling, the 
absolute square of this latter gives the photon number 
per atom). Therefore, in the rest of this paper the sys- 
tem parameters will occur only in the form of the above 
combinations. However, we refrain from introducing a new 
notation for the scaled parameters in order to signify their 
relation to experimental parameters. 

Further simplification for the numerical method can 
be obtained by making use of the periodicity of the op- 
tical potential. We can consider periodic wave functions 
4'{x), and solve Eqs. ^ in the interval [0, A] using peri- 
odic boundary condition. In this way, we discard the dy- 
namics within a "Bloch band" . We are interested in effects 
due to the condensate-cavity field interaction which does 
not couple states of different quasi momenta. Thence it is 
independent of the length of the interval and we choose 
the shortest one to reduce the computational effort. By 
contrast, the coUisional interaction depends on the actual 
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atom density, hence, for a fixed atom number A^, the den- 
sity is artificially enhanced by folding the space into a 
single A interval. Therefore, the collision parameter has to 
be modified such that Ng^X correspond to the collisions 
at the actual (experimental) atom density in the cavity. 



3 Steady state 

In this section, we intend to study the steady state of the 
compound condensate-cavity system. It is described by 
the field amplitude a^, and the condensate wave function 
il}{x^t) = ipQ{x)e~^^* , where hfj, is the chemical potential. 
Following from Eq. (jlj , they obey the system of equations 



N{^o\vt{x)\i^o) 



Ac-N{^o\U{x)\%)+ti^ 



(2a) 



-t- \ao\^U{x) + 2Re{ao}r]t{x) 



+ Ng,\Mx)\'' Uo(x) = fiMx) . (2b) 



\ 2hm 



The solution for iPq{x), ii and ao can be determined in a 
self-consistent manner. We use a variant of the imaginary 
time propagation method. Our method consists of prop- 
agating the condensate wave function in imaginary time 
T = it according to Eq. ( lb I complemented by the adi- 



abatic elimination of the cavity field dynamics. The field 
amplitude is expressed in terms of the instantaneous wave 
function ^'{x), similarly to the form of Eq. (2a I, and in- 
serted into Eq. ( lb I in each time step. Starting from a 
guessed initial function, the wave function '0(x) decays 
with time r. It is just the solution ipo that undergoes the 
slowest decay, because it has the smallest energy S/x. Com- 
ponents from the excited states with higher energy decay 
faster. Therefore, by renormalizing ^'{x) in each step to 
have a norm 1, the iterations converge to the solution 
ipoi^x). This solution is a self-consistent ground state of 
the condensate in the optical potential corresponding to 
ao- 



3.1 Brief summary on cavity-induced self-organization 

The mean-field model of the condensate is closely related 
to the mean-field model of a thermal gas of atoms pre- 
sented in P^EU] . Therefore we can expect that the phase 
transition appears as an abrupt change of the self-con- 
sistent steady state at a critical pumping strength t]. The 
physical process is analogous, however, at zero tempera- 
ture the kinetic energy and the atom-atom collisions are 
apt to stabilize the phase with larger symmetry. In the 
following, we briefly summarize the basic elements of self- 
organization. 

In order to gain some physical insight, it is appropriate 
to adiabatically eliminate the cavity field. The resulting 
self-consistent potential in the GPE ( lb I is 



which is the sum of a A and a A/2 periodic potential. The 
coefficients are 

Ui ^ 2 (cos kx) NIa [Ac - NUo{cos^ kx)] , (4a) 

U2 = {coskxf N'^IqUo . (4b) 

The nonlinearity of the system is introduced by the de- 
pendence of the coefficients on the wave function via the 
specific mean values 



= (^/lol cos kx \ipo) 



(5) 



= 1 and ao = is a 
(cos kx) — for uni- 



which can be considered as an order parameter, and 

B= {'^o\cosHx\^o} , (6) 

which is called the bunching parameter. The order param- 
eter describes the A-periodic spatial order of tjjo- O = 
for a uniform distribution, and = ±1 for the atoms be- 
ing localized around the 'even' {kx — 2mr), or 'odd' field 
antinodes {kx = {2n + 1)tt), respectively. The bunching 
parameter refiects the degree of localization of the atoms 
into the potential wells of the optical lattice. Finally, Iq, 
scaling the depth of the potential, represents the max- 
imum number of photons an atom can scatter into the 
cavity 

° [Ac - NUoBf - 

It can be easily seen that ■00 (a^) 
trivial solution of Eqs. Because 
form atomic distribution, the cavity field driving term van- 
ishes and ao = 0. There is no potential then, Ui = U2 = 0, 
thus the uniform distribution ^o{x) — 1 remains a good 
solution. However, it may not be a stable self-consistent 
ground state of the condensate, because fiuctuations of 
the wave function can be amplified for adequately chosen 
parameters. 

Self-organization is based on the effect of the A peri- 
odic Ui potential term which involves a positive feedback 
mechanism. Let us set the detuning Ac < —N[Uo[ so that 
the sign of Ui becomes the opposite of the sign of 0, cf. 
Eq. ( 4a I . Consequently, if some A-periodic fiuctuation of 
the condensate yields > {i.e., more atoms happen to 
be near the even sites than near the odd ones), it produces 
a potential Ui coskx (with Ui < 0) which has minima at 
the even sites, thus attracting even more atoms there. Si- 
multaneously, the A-periodic optical lattice of condensate 
atoms fulfills the Bragg-condition for constructive interfer- 
ence, {r]t{x)) ^ 0, so the atoms can scatter pump photons 
into the cavity. The positive feedback requires then just 
that the parameters Ui and have opposite signs. This 
can be ensured by choosing the cavity detuning 



Ar 



NUo 



(8) 



V{x) = Ui cos kx + U2 cos^ kx , 



(3) 



which is a sufficient but not necessary condition (because 
the bunching parameter is less than 1, the field could be 
tuned closer to resonance). 

The A/2 periodic U2 potential term does not discrimi- 
nate between the even and odd sites. Moreover, it is pro- 
portional to the square of the order parameter, thus it 
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Fig. 1. The order parameter O, plotted as a function of the 
transverse pump amplitude V^-q, exhibits the self-organiza- 
tion of the BEG into a A-periodic optical lattice. Parameters: 
Ng^ = lOhujRX, NUo = -WOujB., = -iQQujR, k = 200tJfl. 



plays no role in the onset of self-organization. The conse- 
quences of this potential term will be discussed later. 

The runaway solution is counteracted by the kinetic 
energy and the collisions, both trying to maximally spread 
the atomic wave function. There is a critical pump power 
above which the self-organization can occur. 



3.2 Results from the Gross-Pitaeskii equation 

The relevant parameters are measured in units of the recoil 
frequency ojn = hk'^ /(2m) and in units of the recoil energy 
hojn, respectively. The natural length scale is the optical 
wavelength A of the cavity field. We set the strength of 
the s-wave interaction Ngc = 10 uijiX. The atom-cavity 
coupling is characterized by NUq — 100 ujjj, which is on 
the order of the cavity decay rate, k = 200 cj^,. We set 
the value of the cavity detuning, according to Eq. ([s]), 
Ac = -300 Wfl. 

The numerical solution for the self-consistent ground 
state of the condensate-cavity system confirms the quali- 
tative picture presented above. In Fig. [T] the appearance 
of the self-organized lattice is manifested by the variation 
of the order parameter against the pumping strength 

Above a critical pump amplitude {VNr]c ~ 65.6 w^), 
the homogeneous condensate, our initial guess for the con- 
densate wave function, self-organizes into a lattice of pe- 
riod A. At the outset of the self-organization, a sponta- 
neous symmetry breaking occurs: the condensate occupies 
either the kx = site or the x = ir site. Further increasing 
the pump strength -v/TV 77, the atoms get more and more 
localized around the chosen site, that is indicated by the 
growth of the order parameter 0. The localization of the 
condensate in the self-organized phase is showed in Fig. |2] 
for two specific transverse pump amplitudes ^/Nrj =100 
(thick lines) and 300 lur (thin lines) . 

This localization behavior can remind the reader of a 
superfluid-Mott (SM) phase transition (see e.g. Ref. [H] 
for a review) . An SM transition happens in a deep enough 
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Fig. 2. Two typical atomic position distibutions |i/>o(a;)P (solid 
lines}and the corresponding adiabatic optical potentials from 
Eq. (|3| (dashed lines) for pump strengths V~N rj = WOluh (thick 
linesjand \/]V 77 — SOOlur (thin lines) . Other parameters are the 
same as in Fig. [l] 



optical lattice, where the competition between kinetic en- 
ergy and on-site interaction results in distinct ground- 
state configurations. If the on-site interaction dominates, 
the multiple occupation of a single site is energetically 
disfavored and the ground state corresponds to localized 
atoms evenly distributed over the lattice points (Mott 
phase). Contrarily, when the kinetic energy dominates, 
the zero quasi momentum state (completely delocalized) is 
macroscopically occupied (superfluid phase). This shows 
similarity to the self-organization discussed above, with 
the difference that the competition now is between the 
kinetic energy and the energy of the cavity field. This dif- 
ference, however, leads to some other distinctions between 
the two transitions. In the case of the self-organization in 
an cavity the optical potential builds up for r] > rjc and is 
shallow in the beginning. Therefore the macroscopic wave 
function, although it breaks the continuous translational 
invariance, is not strongly localized (at least for rj « rjc). 
Even if the optical potential is deep (77 ^ ijc) the com- 
pressibility of the system, dN/dfi, is nonzero, which is an- 
other difference compared to the SM transition. The role 
of the interparticle interaction strength is also opposed 
in the self-organization transition, since it spreads out the 
macroscopic wave function and counteracts localization. It 
is also important to note, that even in the deeply localized 
limit of the self-organized phase, the number of bosons in 
a single site can be large (depending on the parameters of 
the actual realization) and the Gross-Pitaevskii approxi- 
mation applies. For example, in the experiments made by 
the Esslinger group [7 , the number of condensed atoms 
were approximately iV = 4 • lO'"^, while the size of the con- 
densate were about 30A. Therefore there can be about 10^ 
atoms per lattice site in a possible experimental realiza- 
tion, which is a huge number compared to those found in 
a Mott phase. 

The critical pump amplitude V Nrjc can be determined 
analytically by the stability analysis of the trivial solution 
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of Eqs. being iJo{x) = 1, ao — 0, fiQ — Ngc- Only the 
Fourier component cos kx of a noisy perturbation of the 
wave function can produce a non-vanishing cavity field 
ao other than zero, in accordance with the mean value 
{ijj\cos kx\4>) in the numerator of Eq. (2al. Therefore, 
starting from the wave function ip{x) = 1-1- ecos kx, with 
e 1, we carry out one iteration step of the imaginary 
time propagation method: 



All; 



Ngc — ecos kxx 

2Ac - NUo 



LUR + Nt]' 



{Ac - + 



+ iNgA . (9) 



The component V'o = 1 decays with the rate /iq = Ng^ 
while the decay of the perturbation cos kx depends on the 
pumping strength "q. In order to preserve the stability of 
the homogeneous ground state, the perturbation should 
decay faster than -^o- The condition that the coefficient of 
e cos kx is equal to Ng^ leads to the critical pump ampli- 
tude 



Inte restingly, the critical pump amplitude given by Eq. 
( 10 ) is formally analogous to the one obtained in the case 
of a thermal classical gas [ID] . The temperature is substi- 
tuted by the kinetic energy and the s-wave collision of the 
atoms, formally kBT/h <-> ujji + 2Ngc. 



4 Collective excitations 

Let us now calculate the excitation spectrum of the cou- 
pled condensate-cavity system as the Hnear response of 
the self-consistent steady-state. Two limiting cases can be 
relatively easily understood: (i) for 77 = there is no field 
in the cavity and one gets back the excitation spectrum 
of a homogeneous Bose-gas; (ii) for 77 —^ 00, deeply in 
the self-organized phase, where the optical potential can 
be approximated as a parabola, one expects to obtain the 
excitations of a BEG in harmonic trap potential. This sim- 
plification, however, does not perfectly apply since the ex- 
citations are not only collective oscillations of the atom 
cloud, but they are polaritons involving the fluctuation of 
the field amplitude around its steady state. In the follow- 
ing, we will consider the full transition range, including 
the critical point, between these limiting cases. 

We need to consider the deviations from the stationary 
state (V'o and ao): 

a{t) = ao + Sa{t) , (Ha) 

il;{x,t) = e-'^^'i^x) + S^x,t)] . (lib) 

Inserting the ansatz into Eqs. ([ij and linearizing in Sip 
and Sa, one gets 

iS& ^ ASa + NaoiiM U{x) IStp) + {dip\ U{x) \iPq)] 

+ 7V[(Vo| Vt{x) \S4') + {SiP\ mix) \i>o)] , (12a) 



idijj = {Ho + Ngc\Mx)\^} H + Ngcijl{x)5r 
+'4!Q{x)U{x){aQ5a* + agJa) 
+il}o{x)rjt{x){5a + 5a*) , 



where 



and 



A = - Ac + N {MU{x)\^;o) ~ 



(12b) 



(13a) 



Ho 



2hm 



Ngc\Mx)\'' - ^^ 

+ \ao\^U{x) + 2Re{ao}r7t(a;) • (13b) 



Because the linearized time evolution couples Sip and Sa 
to their complex conjugates, we search the solution in the 
form 



5a{t) 



t X * 



Sip{x,t) = e-"^'Sip+{x) + e'^^'^SipHx) 



(14a) 
(14b) 



where lu = v — is a, complex parameter of the oscilla- 
tion s tand ing for frequency v and damping rate 7. Equa- 
tions ( 12 ) have to be obeyed separately for the e~'"* and 



terms, which leads to the linear eigenvalue equation: 



(15) 



where M is a non-Hermitian matrix being determined by 
Eqs. ( 12 1. Collective excitations of the system are the solu- 





/ Sa+ \ 










= M 




UJ 


SiP+{x) 


H+{x) 




\SiP_ix)J 




W-{x)J 



tions of this eigenvalue problem. It simplifies if we choose 
ipoix) real, and write M in the basis of the symmetric and 
antisymmetric combinations. 



Saa — Sa^ — Sa^ , 
Sag — Sa^ + Sa- , 
S!{x)^SiP+{x) + Si,-(x) , 
Sg(x) — Sip-{x) — Sip^{x) . 

Transformed into this basis, the matrix M is 



(16) 



■ 


-Ho 

(17) 

The argument x has been omitted from the functions ipo, 
U and rjt for brevity. The integral operators X and Y, 
coupling the condensate excitations into the field, read 



~iK ReA 2A^(ReaoX + Y) 

ReA -IK 2iNlmaoX 



2iiPoUlmao -2iPo{UReao + Vt) -Ho - 2Ngcipl 



Xe(x) = / dxil;o{x)U{x)^{x), 



Y^(2:) = / dxipo{x)'qt{x)£,{x). 



(18a) 



(18b) 



Being a non-Hermitian matrix, M has different right and 
left eigenvectors corresponding to the same eigenvalue w. 
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The matrix has a special symmetry: if lo is an eigen- 
value with the right eigenvector 

{5aa,5as,5f, Sg) , 

then —uj* is also an eigenvalue with the right eigenvector 

{-6al,Sal,Sf*,-5g*) . 

Thus the eigenvalues come in pairs, having the same imag- 
inary parts but the real parts are of the opposite sign in 
a pair. This grouping of the eigenvectors makes sense for 
eigenvalues with non- vanishing real part. 

The imaginary part describes damping which arises 
from the nonadiabaticity of the cavity field dynamics. The 
linear perturbation calculus of the excitations which we 
adopted here takes into account that the field follows the 
changes of the BEC wave function with a delay in the or- 
der of Depending on the specific choice of the param- 
eters, this can yield damping or heating, which is known 
as cavity cooling and has been extensively studied in the 
past decade. Damping of BEC excitations in optical cavi- 
ties, which has been first studied in Refs. [T5(l22j. is an in- 
teresting opportunity which motivates the realization and 
the study of the coupled cavity-condensate system. 

4.1 Spectrum analytically below threshold 

Below threshold, the condensate wave function is constant 
and linearization around this simple solution leads to an- 
alytical results. The stability of the uniform distribution 
can be evaluated from the spectrum, thus the critical point 
can be determined in this way, independently of the previ- 
ous calculation in Sec. 



3.2 



Moreover, we obtain a detailed 
description of a restricted part of the spectrum which con- 
tains polariton excitations. 

Assuming homogeneous atomic distribution and, cor- 
respondingly, vanishing cavity field (i/'o = 1; Q^o — 0); 
the only non-trivial coupling term in the matrix M, see 
Eq. (17 1, is the Y operator defined in Eq. (18d). Since 
77t(x) cx cos /ex, only the Fourier component cos kx cou- 
ples to the cavity amplitude. In return, the fourth line of 
the matrix M shows that the field fiuctuations excite just 
this rjt^x) oc cos fca; condensate perturbation. Thus, this 
subspace is closed below threshold. All the other conden- 
sate excitations decouple from the field and remain simply 
the higher Fourier components with frequencies 



^n?ujR{v?ujR + 2Ngc) , {n > 1) 



(19) 



which is identical to the excitations of a homogeneous con- 
densate in a box [53]. It is sufficient to diagonalize the 
matrix M in the restricted subspace. 
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It has the fourth order characteristic equation 
(A^ - f2f) [{iK + Xf - Si] ~ 2Nri'LOR5c = 



(20) 



0, (21) 
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Fig. 3. Eigenvalues of the lowest condensate excitation of a 
homogeneous BEC from Eq. (211 as a function of the pump 
strength V^-q. As shown in the inset, the real parts (solid red) 
vanishes slightly below the critical point which is reached when 
the upper branch of the imaginary parts (dashed green) crosses 
zero. In the main figure the imaginary parts are magnified by 
a factor of 50. Parameters are the same as in Fig. [l] 



where Sr 



Jc = —^c + ^NUq, and f2i is the first excitation 
energy in the Bogoliubov spectrum ( 19 1 of a BEC in a box. 
One can check that, for = 0, the last term vanishes and 
Ai,2 = ±f2i for the condensate excitation, A34 = zLSc — in 
for the cavity mode. For non-zero pumping, 77 0, the 
condensate-like and the field-like excitations mix. When 



hNUn 



nl « 



which case we consider in this paper, the 
frequencies corresponding to the excitations of the free 
field and to the free condensate are well separated. There- 
fore, the mixing ratio is small so that the polariton modes 
can be attributed to dominantly condensate or field exci- 
tations. 

At the onset of self-organization the uniform ground 
state "00 changes, which, in the present approach, cor- 
responds to the appearance of a zero eigenvalue in the 
spectrum. In Fig. [3] we plot the numerical solutions of 
Eq. (21 1 for the lowest, dominantly condensate- type ex- 



citation. The real parts (solid red) tend to zero as in- 
creasing the pump strength. Oppositely, the absolute value 
of the imaginary parts of the eigenvalues (dashed green) 
increases with the pump strength. The behaviour of the 
eigenvalues near the critical point y/Nrjc — 65.612 0^7? is 
magnifiedjn the inset of Fig.|3] When the real parts reach 
zero at V Nrj^, , the initially identical imaginary parts split 
up and the upper branch crosses zero. This crossing is the 
critical point, here the "00 steady-state becomes dynami- 
cally unstable. By expressing 77 from Eq. (21 1 for A = 0, the 



same critical transverse pump amplitude r]c is obtained as 
that in Eq. pB. 



The emergence of a positive and a negative imaginary 
part just at the vanishing of the real parts of an excita- 
tion is a typical feature of a condensate's instability. In our 
case, however, coupling to the cavity field yields a negative 
imaginary part where the real part becomes zero. Thus, 
there is a narrow range, as shown by the inset of Fig. [3] 
where the real parts remain zero, but both imaginary parts 
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are negative. This unfamiliar course of the dynamical in- 
stability is a signature of cavity cooling, thus it is so for 
finite cavity decay rate k. 



In the limit Q'f < + 5}., 
quency can be approximated by 



the first excitation fre- 



ReAi — f2i\ 



(22a) 



and the imaginary part is quadratic in the pump strength, 



(22b) 



These approximate expressions fit well the curve shown 
in Fig. [3] The s-wave interaction between atoms increases 
the decay rate of this particular excitation. The imaginary 
part vanishes, of course, for k ^ 0. 

The width of the narrow range, where the V'o = 1 uni- 
form steady state is stabilized by that the cavity cooling 
damps out zero-energy excitations, is obtained: 



2 2 2 



K f2l 



(23) 



Similarly to the decay rate ImAi, the range expands on in- 
creasing the coUisional parameter Ngc- Finally, note that 
Eq. ( 23 1 is considered a small parameter in the estimation 
of Ai in Eqs. (p2|. 
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4.2 Full spectrum 

Above threshold, the collective condensate-cavity excita- 
tions belong to the self-consistent ground state given by 
Eqs. PI which can be calculated only numerically. There- 
fore, the solution of the eigenvalue problem in Eq. ( 15 1 



is performed also numerically using LAPACK. The wave 
function is defined on a spatial grid of 200 points in the 
interval of one wavelength. 

In Fig.Qa), we show the first three excitation frequen- 
cies and damping rates (71/2 is plotted) of the conden- 
sate around the real ground state, and the frequency and 
damping rate of a higher excited state corresponding to 
the cavity mode (double dashed lines, and 7y/4000 

are plotted). We depict only the positive frequencies, tak- 
ing into account the symmetry of the eigenvalues of the 
matrix ([l7]). 

Below threshold the wave function ipoi^x) is constant, 
only the first excitation cos kx couples to the cavity field, 
and tends to zero on increasing the pumping, as it has 
been discussed before. All the other excitations are inde- 
pendent of the field, and have constant energies giving 
back the Bogoliubov excitation spectrum of a condensate 
in a box, see Eq. (19). Two orthogonal excitations have a 



fixed number of nodes and are degenerate in this regime. 
One of the modes in each pair is orthogonal to the cavity 
mode function cosfca;. 

At the onset of self-organization the eigenvalues dra- 
matically change. At the critical point the excitation en- 
ergies drop slightly below the Bogoliubov energies given 



Fig. 4. Frequencies (a) and decay rates (b) of the six lowest 
collective condensate excitations and the ones corresponding to 
the cavity field as a function of the transverse pump strength. 
Decay rates of the lower branches of each pair are 0. For the 
cavity-dominated mode the frequency and the decay rate is 
divided by 5 and 4000, respectively, and 71 is divided by 2 in 
(b). Parameters are the same as in Fig. [l] 



in Eq. ( 19 ). This dip is related to the collisions and disap- 
pears for — > 0. Let us mention that the numerical cal- 
culation of the stationary solution of the Gross-Pitaeskii 
equation becomes inaccurate at the critical point. The 
convergence of the iterative solution slows down, which 
is an inherent consequence of the criticality obtained at 
the degeneracy of the ground and first excited state en- 
ergies. Nevertheless we checked in some points that the 
dip is still there if very high accuracy is demanded in the 
iteration process. 

At threshold, the degeneracy in the excitation pairs is 
lifted. The lower branch remains orthogonal to the cavity 
mode and decouples from the field fluctuations. The up- 
per branches correspond to polariton excitations mixing 
condensate and field fluctuations. 

Far above threshold, ^/Nr] 00, the excitation fre- 
quencies increase linearly with ^/N t] and are uniformly 
spaced. This can be understood if we approximate the 
deep optical trap by a harmonic potential in the vicinity 
of the antinodes. The adiabatic potential in Eq. ([3]), trans- 
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formed by using cos fcx w 1 — (fcx)^/2, gets a characteris- 
tic harmonic frequency proportional to the square root of 
the intracavity intensity. Thus, from Eq. ([t]), the trap fre- 
quency is linearly proportional to ^/Nrj. The spectrum is 
composed of the integer multiples of this frequency. This 
is obvious in the coUisionless case, Qc = 0, where the exci- 
tation frequencies are the same as for a single particle in a 
harmonic potential. The other extreme case, the Thomas- 
Fermi regime which we are closer to with the parameters of 
Fig. |4] can also be treated. Here, for a three-dimensional 
harmonic trap the excitation frequencies are given by a 
linear combination with integer coefficients of the three vi- 
brational frequencies |24l . This can be contracted to one 
dimensional motion by assuming very large frequencies 
in the transverse directions, then the low excitations are 
obtained as integer multiples of the longitudinal trap fre- 
quency. 

In Fig. [4jD, the decay rate of the excitations corre- 
sponding to the upper branches in Fig. 4|a) are shown. 
The lower branch of each pair has zero clamping because 
it is orthogonal to the mode function of the cavity field. 
The lowest excitation decays with a rate 71 which exhibits 
a dip at threshold. In principle, it should drop down to 
zero, as shown in Fig. [3] representing the exact result be- 
low threshold. The present result, based on the numerical 
calculation of the ground state, is not accurate enough 
in the vicinity of the threshold to resolve the vanishing of 
7i . However, somewhat above threshold the numerical ap- 
proach becomes very accurate. It shows then that 71 rises 
back to the value it had below threshold and then increases 
further with increasing y/Nrj. Not too far above threshold, 
the weakly localized ground state enables the coupling of 
higher condensate excitations to the cavity field. These 
modes become damped, but their damping rate vanishes 
in the y/Nrj 00 limit. There, as the ground state tends 
to a tightly localized one in a harmonic trap, only the 
second excitation is coupled to the field [15] . 

The frequency and the decay rate of the cavity field 
are weakly perturbed by the condensate. The frequency 
of the cavity mode is expected to be i^/ = Ac — NUqB, 
that depends on the collective coupling NUq and, through 
the bunching parameter, on the localization of the ground 
state ipQ. Decrease in the field mode frequency is accom- 
panied by an increase of the stationary photon number Iq 
in the cavity. 



5 Strong collective coupling regime 

The interpretation of self-organization in terms of adia- 
batic optical potentials in Sec. |3.1| reveals that the effect 
relies on the discrimination between the sites kx = and 
kx = TT (even and odd sites). The symmetry breaking is 
attributed entirely to the A periodic potential. The other, 
A/2 periodic cos^ kx potential could have been disregarded 
in this respect. Above threshold, however, this latter po- 
tential can have a significant role. If the condensate local- 
izes around, say, kx = 0, the cos kx potential has max- 
ima a,t kx = n and atoms are repulsed from this region. 



which assists the self-ordering process. The cos^ kx poten- 
tial term has always minima both at fca; = and kx = tt, 
thus it counteracts the repulsion. For certain parameters 
it may occur that a secondary minimum appears and the 
adiabatic potentials form an asymmetric double well po- 
tential. So called 'defect' atoms scattering photons with 
'wrong' phase can be trapped in the shallower traps. 




50 100 150 200 250 300 
\/N 7] [in units of c^tj] 



Fig. 5. Full phase diagram of the system for collision param- 
eters Ngc — and 10 Xlur. For Nqc = the dashed blue 
curve represents the phase boundary between homogeneous 
distribution and self-organized lattice. In the region between 
the dashed blue and the solid red line (these curves coalesce 
asymptotically for large NUo) there is a self-organized lattice 
without defects. The solid red curve gives the critical NUo be- 
low which the secondary minima occur. Right shifted, the same 
phase boundary lines are drawn by dotted green and dashed- 
dotted brown for Ngc — WXuir. Parameters are k — 200ujr, 
Ac = -2k. 



In the organized pattern, the emergence of a secondary 
potential minima at the complementary sites depends on 
the ratio of the two potential terms in Eq. ([s]) , namely on 
Ui/U2- For perfect localization (rj — > 00), the condition 
for the possibility of stable defects is 



N\Uo\ > K 



(24) 



which is precisely the condition of strong collective cou- 
pling in cavity QED. For finite pumping strengths {Nrj < 
00), the ratio of the two potential terms is dependent 
on the order parameter 0, thus on the localization of 
the condensate wave function. Therefore, the above con- 
dition Eq. p4| ) should change with the pump strength 
Nrj and with the s-wave collision parameter gc- Because 
U2/U1 oc < 1, the larger the order parameter is, 
the smaller A^|J7o| is needed for getting defects. The crit- 
ical lines are plotted in Fig. |5] below which the secondary 
potential minima appear as a function of the pumping 
strength for two different values of the s-wave collision pa- 
rameter Qc = and Ngc = lOXcufj. The cavity detuning is 
fixed to Ac — —2k. When collisions are not negligible, one 
needs a larger light shift iV|[/o| to reach the regime where 
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defects are expected, just because the collisions reduce the 
order parameter O. The vertical asymptote of the curves 
is at the threshold of the self-organization, while the hor- 
izontal one is at the same value defined by Eq. ( 24 1 that 
describes the \/Nrj oo behaviour for perfect localiza- 
tion. 
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Fig. 6. Frequencies (a) and decay rates (b) of the six lowest 
collective condensate excitations in the regime, where the sec- 
ondary potential minima are present as a function of the trans- 
verse pump strength. Decay rate of the first excitation 71 (solid 
red) is divided by 40. Parameters are gc — 0, NUo = —IOOOcjr, 
K = 200 and Ac = -1200 wh. 



Now let us turn to the effect of the secondary poten- 
tial minima on the collective excitation spectrum. In the 
regime where the energy difference between the two types 
of potential minima at the complementary sites is in the 
order of their depths, the lowest excitation energies in the 
deeper wells are comparable to those in the secondary 
ones. There will be collective excitations that are com- 
posed of the combination of excitations localized in the 
two different wells. 

For simplicity, we calculate the spectrum in the coUi- 
sionless case for = 0- In Fig. [6] we plot the collective 
excitation frequencies (a) and the decay rates (b) as a 
function of the pumping strength i] for k = 200w/j and 




Nrj/ujR 



Fig. 7. The 5tl)+{x) components of the eigenvectors of the 
collective excitations with the lowest (a) and the highest fre- 
quencies (b) in Fig. are plotted for chosen values of rj. In 
these special cases the other components are zero. 



NUq = — IOOOlur. Starting from the Bogoliubov spectrum 
of the homogeneous condensate, there appear two types 
of excitation modes in the self-organized phase. The exci- 
tations that are localized in the deeper wells are the same 
as the ones presented in Fig. |4] Notice the lack of the dip 
at the critical point, which is because of the choice gc = 0- 
The other type, represented by the two excitation modes 
bending up in Fig. contain the excitations of the 'de- 
fect' well. For perfect localization, both the lower and the 
higher potential minima, hence also their difference, are 
proportional to N\r]\'^. Exciting the condensate into the 
defect well costs the energy difference, this is the reason 
of the quadratic behavior as a function of the pumping 
strength y/N rj. 

In order to exemplify the effect of the secondary po- 
tential minima in the spectrum, we plot the wavefunctions 
of two condensate excitations at increasing values of the 
pump strength y/N rj in Fig. [rj Namely, we chose the ex- 
citations with the lowest and the highest frequencies in 
the right of Fig. [6^. Both of them are decoupled from the 
cavity field so that the wavefunction is real and can be in- 
terpreted as a condensate wavefunction in position space. 
The lowest one (a) is proportional to sin kx in the uniform 
phase at \/Nri — 15, having two nodes at a; = and A/2. 
On increasing the pump strength y/Nrj, the condensate 
wave function gets more and more localized, hence this ex- 
citation contracts into the primary potential well, which 



10 



D. Nagy, G. Szirmai, P. Domokos: Self-organization of a Bose-Einstein condensate in an optical cavity 



is now centered at x = A/2 in the self-organized phase. 
At y/Nrj — 120, the strongly localized BEC feels just the 
harmonic term of the optical potential. Therefore, at this 
end, we get for the wavefunction of the excitation some- 
thing close to the first excitation of a harmonic oscillator. 
The other selected excitation is the highest one, bending 
up in Fig. |6|i. It has six nodes in the homogeneous phase, 
that corresponds to the third excitation of the Bogoliubov 
spectrum. It runs side by side with its orthogonal pair up 
to y/Nrj « 60, where they split up. The lower branch tends 
to the fourth harmonic oscillator excitation in the primary 
well for y/Nrj oo, however the upper branch plotted in 
Fig. [TJd becomes the first excitation in the secondary well 
centered at a; = 0. 



6 Quantum depletion 

A crucial issue on the steady-state solution of the coupled 
Gross-Pitevskii equations (jlj (besides the linear stabil- 
ity analysis of the previous sections) is the magnitude of 
the quantum depletion which is a small parameter in our 
model theory. In our situation depletion is caused by two 
physical processes. 1) The atom-atom interaction kicks out 
atoms from the condensate, such as in a conventional in- 
teracting Bose gas. 2) The atom-light interaction mani- 
festing itself in the polariton like resonances prevents the 
simple product wavefunction of the Gross-Pit aevskii the- 
ory to be a true eigenstate of the coupled system. To syn- 
thesize the quantum depletion caused by the polariton na- 
ture of some of the excited states, we perform its calcula- 
tion with the interatomic collision set to zero {gc = 0) . We 
restrict this study to the case of a lossless cavity (k = 0), 
so that the system be conservative and the number of non- 
condensed atoms {N') can be calculated in the usual way, 
based on the eigenvectors of the Bogoliubov matrix M 
[23]. 

Without the detailed calculation, we present the re- 
sult for the quantum depletion in Fig. [8] The sharp peak 
in the quantum depletion indicates the critical point. It is 
not surprising that the depletion increases up at the criti- 
cal point where the first excited state becomes degenerate 
with the ground state. In fact, the number of atoms in the 
first excited state diverges as the pumping t] tends to 
the critical point from below as 



8Ai 



1 



8 \SN{j^,-r,y 



(25) 



where we used Eq. (22a I to approximate the first eigen- 
value Ai near the threshold. The mean field theory and 
the perturbation approach is valid if the depletion rate 
N' /N is much smaller than 1. With reasonable experi- 
mental values of TV, i.e. ^ 1, this condition is obeyed 
with the exception of a very narrow range around the crit- 
ical point. This justifies the use of the mean field model 
in almost the full range of Fig. |8] Moreover, by means 
of the divergence, the present model is suitable to locate 
precisely the critical point. More details on the quantum 
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Fig. 8. Quantum depletion due to the coupling to the cav- 
ity field as a function of the pumping strength in a semi- 
logarithmic plot. A'^' is the number of atoms distributed in the 
excited states, for fractional depletion the plotted number has 
to be divided by the actual number of atoms TV. Collisions 
are discarded, Qc = 0. Other parameters: k = 0, Ac = —300, 
NUo = — 10 (solid red curve) and NUo — —100 (dashed green). 



depletion and its treatment in an open system with k =/= 
will be the subject of a forthcoming paper. 



7 Conclusion 

Quantum many-body systems formed by neutral atoms 
weakly interacting through the electromagnetic field is of 
general interest. One possibility consists in using a dipolar 
gas [15] • Although, the dipole-dipole interaction is rela- 
tively weak and of short range, patterns, collective modes 
[26 and instabilities [27] arise in a dipolar condensate. The 
other possibility is that the atoms with tiny induced dipole 
moment are placed inside a cavity and couple strongly to 
one given mode of the resonator. Then the interaction 
has a long range, all atoms are coupled to all the other 
within the cavity. In this system, for example, the real- 
ization of the Dicke-model quantum phase transition is 
proposed and seems feasible [2511^150] by means of using 
Raman- type transitions in multilevel atoms. Another ex- 
ample, which concerns the motional degrees of freedom of 
particles, is that a thermal gas of laser-driven atoms pro- 
duces a classical phase transition called 'self-organization' 
in a cavity [I]. In this paper we addressed the quantum 
version of this latter problem and gave a detailed account 
of the steady-state and the dynamics when the atoms form 
a Bose-Einstein condensate in the cavity. 

Within the framework of a Gross-Pitaevskii like mean- 
field model, we showed that the steady-state of the driven 
condensate is either the homogeneous distribution or a A- 
periodic ordered pattern, and the two regimes are well sep- 
arated by a critical point. That is, the quantum analogue 
of the classical self-organization phase transition exists for 
Bose-condensed atoms. The critical point, corresponding 
to a threshold pump power has been determined analyti- 
cally from the Gross-Pitaevskii equation, and also analyt- 
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ically from the spcctrmn of the excited states. We showed 
that the spectrum of the collective excitations becomes 
more intricate in the strong collective coupling regime of 
cavity QED, where the ordered phase is described by an 
asymmetric double well potential and defect atoms can 
appear. We showed that the quantum depletion diverges 
in the critical point. Otherwise, outside a narrow range 
around it, the coupled condensate-cavity system is stable 
and the mean-field approach is well justified. 

We acknowledge funding from the National Scientific 
Fund of Hungary (NF68736, T046129, and T049234). 
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